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The inviscid Burgers equation with random and spatially smooth forcing is considered in the limit 
when the size of the system tends to infinity. For the one-dimensional problem, it is shown both 
theoretically and numerically that many of the features of the space-periodic case carry over to 
infinite domains as intermediate time asymptotics. In particular, for large time T we introduce the 
concept of T-giobal shocks replacing the notion of main shock which was considered earlier in the 
periodic case (1997, E et al., Phys. Rev. Lett. 78, 1904). In the case of spatially extended systems 
these objects are no anymore global. They can be defined only for a given time scale and their 
spatial density behaves as p(T) ~ y- 2 / 3 f or j ar g e j\ The probability density function p(A) of the 
age A of shocks behaves asymptotically as A~ 5 ^ 3 . We also suggest a simple statistical model for 
the dynamics and interaction of shocks and discuss an analogy with the problem of distribution of 
instability islands for a simple first-order stochastic differential equation. 
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I. INTRODUCTION 

The d-dimensional forced Burgers equation 

d t u+(u ■ V)« = vV 2 u~VF(x,t), u = -\7ip (1.1) 

which describes the dynamics of a stirred, pressure less 
and vorticity-free fluid, has found interesting applications 
in a wide range of non-equilibrium statistical physics 
problems. It arises, for instance, in cosmology where it 
is known as the adhesion model in vehicular traffic 
H or in the study of directed polymers in random media 
[3j. The associated Hamilton-Jacobi equation, satisfied 
by the velocity potential if>, 

d tV J-^\Vvj\ 2 = u^ + F(x,t), (1.2) 

has been frequently studied as a nonlinear model for the 
motion of an interface under deposition : when the forc- 
ing potential F is random, delta-correlated in both space 
and time, Eq. 11. 2|) is the well-known Kardar-Parisi- 
Zhang (KPZ) equation p|. The case with large-scale 
forcing was considered in Refs. 0, |(| as a natural way 
to pump energy in order to maintain a statistical steady 
state. Burgers equation is then a simple model for study- 
ing the influence of well- understood structures (shocks, 
preshocks, etc) on the statistical properties of the flow. 
As it is well known, Eq. I|l.l|) in the limit of vanishing 
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viscosity (v — > 0) displays after a finite time dissipative 
singularities, namely shocks, corresponding to disconti- 
nuities in the velocity field. In the presence of large-scale 
forcing, it was recently stressed for the one-dimensional 
case 0, @ and also for higher dimensions 0, [ljj, that 
the global topological structure of such singularities is 
strongly related to the boundary conditions associated to 
the equation. More precisely, when, for instance, space 
periodicity is assumed, a generic topological shock struc- 
ture can be outlined. It plays an essential role in un- 
derstanding the qualitative features of the statistically 
stationary regime. 

So far, the singular structure of the forced Burgers 
equation was mostly investigated in the case of finite-size 
systems with periodic boundary conditions. It is however 
frequently of physical interest to investigate instances 
where the size of the domain is much larger than the 
scale, so as to examine, for example, the role of Galilean 
invariance |Tlj . The main goal of the present paper is to 
describe the singular structure of Burgers equation 11. 1|) 
in unbounded domains. To achieve this goal we consider 
the case of spatially periodic forcing with a forcing scale 
Lf much smaller than the system size L. More precisely, 
we assume that the force is random and white in time, 
homogeneous, isotropic, smooth and L-periodic in space 
and investigate the limit L/Lf 3> 1, while keeping a con- 
stant injection rate of energy. Hence, for a fixed size of a 
system L we consider a stationary regime corresponding 
to the limit t — ► oo and then study the limit L — > oo with 
the C 2 norm of the forcing growing like L. As we will see, 
most of the concepts associated to infinite time asymp- 
totics introduced for the periodic case can be generalized 
in this limit as intermediate time asymptotics. 
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It is important to mention that the problem we are 
discussing has also another interpretation which has no 
direct relation to Burgers equation. In fact, we are study- 
ing the structure of singularities for variational problems 
associated to generic time-dependent Lagrangians in un- 
bounded domains. From this point of view the random 
forcing is just a natural way to characterize generic time- 
dependence. 

The paper is organized as follows. In Sec. [H] we give 
a brief exposition of the theory in the case of periodic 
forcing. In particular, we introduce the variational ap- 
proach to Burgers turbulence and present results on the 
existence and the uniqueness of the main shock and the 
global minimizer. In Sec. IIIII we introduce the notion of 
T-global shocks and describe their behavior in the time 
asymptotics Lf/u lms <C T <C L/u rms , where w rms is the 
root-mean-square velocity. We also discuss theoretical 
and numerical results for the density of T-global shocks 
and the probability density function (PDF) of the age 
of shocks. In Sec. IIVI we discuss an analogy with sta- 
bility theory for one-dimensional stochastic ordinary dif- 
ferential equation (ODE) and present a simple statistical 
model for the interaction of shocks. Section IV! contains 
concluding remarks and, in particular, an interpretation 
of the results in terms of an inverse cascade. Most of the 
results of this paper are focusing on the one-dimensional 
case. We finish with a brief discussion of possible exten- 
sions to the multi-dimensional case. 



II. PERIODIC FORCING 

In the case of smooth spatially periodic forcing potential 
F(x, *), it was shown in the one-dimensional case @ and 
for higher dimensions that a statistical steady state 
is reached at large times by the solution to the Burgers 
equation. Taking the initial time at -co, the velocity 
field is periodic in space and uniquely determined by the 
realization of the forcing. At an arbitrary time *, this 
unique solution can be expressed in terms of the following 
variational principle [T^.ll3lll3 | involving the Lagrangian 
action 4 : 



Since the force is potential and spatially periodic, the 
mean velocity 



u(x,t) = V x min [.4(7,*)] 

■y:-((t)=x 



4(7,*) = 



i|7(r)| 2 -T( 7 (r),r) 



(2.1) 
dr. (2.2) 



The minimum in 12.1(1 is taken over all absolutely con- 
tinuous curves 7(t) of M. d with r 6 (— oo,*] such that 
7(*) = x. A curve minimizing the action in 1|2.1[) is called 
a one-sided minimizer. It corresponds to a fluid particle 
trajectory which is not absorbed by shocks till time *. 
For all times r < *, such a minimizer is a trajectory of 
the dynamical system corresponding to the Lagrangian 
flow and thus satisfies the Newton equation 



b = 



1 
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[0,L] d 



u(x, t)dx 



(2.4) 



is the first integral of (|1.1|) . When the initial time is 
taken at — oo, so that the statistically stationary regime 
is reached, the mean velocity b is the only information 
remaining from the initial condition. For simplicity, we 
consider in the sequel the case of a vanishing mean veloc- 
ity (b = 0), but all the results remain true for arbitrary b. 

It is easy to show that for Lebesgue almost all x there 
exists a unique one-sided minimizer. The locations where 
there are more than one minimizer correspond to shock 
positions. Those are exactly the positions where the ve- 
locity potential ip(x,t), which is a Lipschitz function, is 
not differentiable, so that the piecewise continuous veloc- 
ity field u(x,t) — —Vip(x,t) has jump discontinuities. 
As * — > — oo, all the one-sided minimizers which origi- 
nated at time * converge backward-in-time to the trajec- 
tory of the unique fluid particle which is never absorbed 
by a shock. This trajectory, denoted 7 gm (-)> is called the 
global minimizer because it is an action-minimizing tra- 
jectory at any time. The statement about uniqueness of 
the global minimizer holds under some natural conditions 
on the forcing potential (see, e.g., Refs. H,E|). Shocks are 
generically born at some finite time and then grow and 
merge. There exist a particular shock structure, called 
the main shock (d = 1) or the topological shock (d > 1) 
which has always existed in the past. This shock can 
be constructed by unwrapping the picture to the whole 
space M d (see Fig. ^for the case d — 1). We then ob- 
tain a lattice of periodic boxes, each of them containing 
a periodic image of the global minimizer and the topo- 
logical shock corresponds to the set of appositions giving 
rise to several minimizers that approach different global 
minimizers belonging to different images of the periodic 
box. 

We focus now on the one-dimensional case with a 
space-periodic forcing potential of period L. In this case, 
the global minimizer is a hyperbolic trajectory of 12. 3|) 
and is thus associated to a smooth unstable manifold 
which is a smooth curve in the position-velocity phase- 
space {x, u). The main shock is defined as the unique po- 
sition giving rise to the left-most and the right-most one- 
sided minimizers approaching the global one backward- 
in-time. By definition, the trajectory of the global mini- 
mizer 7 gm (*) and the main shock trajectory A ms (*) never 
intersect, so that they have to satisfy A ms (t) < 7 gm (*) < 
X ms (t) + L. Hence, the large-scale displacement of the 
two-sided minimizer is the same as for the main shock. 
The Lagrangian flow l|2.3[) is a second-order stochastic 
ODE. Assuming the vanishing of the mean velocity b, 
one may be tempted to think that the displacement of 
a typical trajectory 7(*) scales like the integral of the 
Brownian motion : 



7 ( T ) = -VF( 7 (T),r). 



(2.3) 



(7(*)-7(0)) 2 )«|*| 3 when |*| 



(2.5) 
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FIG. 1: Sketch of the unwrapped picture in space-time for 
d = 1; each minimizer converges backward in time to one of 
the periodic image of the global minimizer. The main shock 
is defined as the location from which emanate two minimizers 
approaching different images of the global minimizer. 



However this scaling does not apply to the minimizers 
because the particularity of such trajectories cannot be 
ignored. Indeed, when investigating the typical displace- 
ment of the global minimizer, we find that 

((7gm(*) ~ 7gm(0)) 2 ) cx |t| when \t\ -> oo. (2.6) 

This result was obtained by numerical study of the be- 
havior at large times of the position of the main shock. 
This gives the large-time typical displacement of the po- 
sition of the main shock which, as we have seen above, 
coincides up to a constant with the displacement of the 
global minimizer. For the numerical investigation, we 




FIG. 2: Average over 500 realizations of the position X(t) of 
the main shock as a function of time for 6 = and a space- 
periodic forcing with period L = Lf. 

chose a forcing which is the sum of independent random 
impulses acting at discrete "kicking" times ^jj . Between 



successive kicks, the solution to the Burgers equation de- 
cays. This method is particularly efficient for obtain- 
ing the solution in the limit of vanishing viscosity, as 
one can use between kicks the fast Legendre transform 
method |l6j| which is directly connected with the varia- 
tional principle (|2.1[) and the Lagrangian picture of the 
flow. The position of the main shock is obtained by a 
Lagrangian method which is based on the analysis of the 
right-most and the left-most positions from which the 
two-sided minimizer is approached. Figure [3 shows the 
diffusive-like behavior of the main shock trajectory for a 
Gaussian forcing confined to the first three Fourier modes 
in the case when the mean velocity b vanishes. 

The diffusive behavior of the global minimizer has 
also a simple theoretical explanation. Indeed, the mean 
square displacement of the global minimizer can be ex- 
pressed through its velocity as 



(7gmW-7gm(0))- 




7gm(s) ds 



(2.7) 



where the velocity of the global minimizer 7 gm (s) is a sta- 
tionary random process. The hyperbolicity of the global 
minimizer implies that the time correlations for this pro- 
cess decay exponentially. Hence, there are just two possi- 
bilities: the mean-square displacement (|2.7II either grows 
linearly with time, which gives 12. 6f) . or it stays bounded 
in a limit t — > oo. The latter behavior can be easily ruled 
out since the trajectory of the global minimizer does fluc- 
tuate. 



III. STATISTICS OF T-SCALES IN LARGE-SIZE 
SYSTEMS AND PDF FOR THE AGE OF 
SHOCKS. 

As we have discussed above, the existence of the main 
shock in the spatially periodic situation follows from a 
simple topological argument. The main shock is unique 
when the system is considered on the torus correspond- 
ing to its period (circle for d — 1). However, if the sys- 
tem is considered on the whole real line (universal cover), 
then there exists, at each moment of time, an L-periodic 
one-dimensional lattice of main shocks. Each of these 
shocks is associated to two minimizers (left-most and 
right-most) which also form a periodic structure. In the 
case of the main shocks the distance between these two 
minimizers tends to L as t — > — oo. For all other shocks 
the corresponding distance tends to 0. In this sense, the 
main shocks are the only global shocks present in the 
system. There is also another way to characterize the 
main shock. Namely, the main shock is the only shock 
which existed forever in the past, that is, the main shock 
is infinitely old, contrary to all other (local) shocks, all 
of them being created at a finite time and hence having 
a finite age. In other words, local shocks can be traced 
backward in time only for a certain finite interval of time. 
The length of this time interval is what we call the age 
A of the shock. 
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When we drop the periodicity condition and consider 
large-size systems, the main shock disappears. But one 
can still introduce the notion of T '-global shocks, that is 
shocks which behave like the main shock when observed 
during a time scale of the order of T. As in the case of pe- 
riodic boundary conditions, there are now also two equiv- 
alent ways to characterize T-global shock. Geometrically, 
one can again consider left-most and right-most mini- 
mizers associated to a given shock. After tracing them 
backward in time for time intervals longer than a cer- 
tain "correlation time" , these two minimizers are getting 
close and start converging to each other exponentially 
fast. T-global shocks are defined as shocks for which this 
"correlation time" is larger than T. This means that for 
backward times i«r, one has d(t) > Lf, where d(t) is 
the distance between minimizers. Equivalently, one can 
say that T-global shocks have a finite age larger than T. 
Below, we study the statistics of the T-global shocks. 

We start with a numerical investigation of the forced 
Burgers equation in a domain of a size L, much larger 
than the forcing length-scale Lf. The force is taken to 
be Gaussian, statistically homogeneous, white-noise in 
time and characterized by its covariance 



(f(x+e,t + T)f(x,t))=B(e)6(T), 



(3.1) 



where B is a smooth, L-periodic function. Lf is the cor- 
relation length, chosen such that B{£) is concentrated at 
scales I < L f <SC L. The injection rate of energy per unit 
length is kept constant (i.e. we choose the space corre- 
lation such that B(0) <x L). Such a hypothesis is essen- 
tial because we found that it prevents any L dependency 

1/2 

of the root-mean-square velocity u rms = (u 2 {x,i)) 
We will now study the behavior of the statistically sta- 
tionary solution for the intermediate time asymptotics 
Lf/u lms <t< L/u lms . For very large but finite L, the 
problem is well-defined and we can suppose that the ini- 
tial condition is taken at — oo, so that a statistical steady 
state is achieved for the solution of the Burgers equation. 

In order to obtain a sketch of the behavior of the 
solution, the limit of infinite aspect ratios L/Lf is in- 
vestigated numerically using a kicked forcing and the 
fast Legendre transform method (see Sec. |nj. More- 
over, we consider here the case of a forcing whose spec- 
trum is concentrated in the wavenumber intervals \k\ € 
[nLf /L, 3ttL / /L] . As we will see in Sec . El the functional 
form of the forcing space correlation B(£) at large scale 
I Lf plays an important role. Numerical observations 
suggest that, at any time in the statistical steady state 
(obtained after sufficiently long integration), the shape 
of the velocity profile is locally similar to the order-unity 
aspect ratio problem, duplicated over independent inter- 
vals of size L f (see Fig. |3(a)| l . More particularly, when 
tracking the trajectories of fluid particles backward in 
time, one can see that the convergence o f the minimizers 
to each other is far from uniform. Figure 3(b) shows that 
the minimizers are forming different branches which are 
converging to each other backward in time, defining in 
space time a tree-like structure. Eventually, the number 




(a) 




FIG. 3: (a) Upper: snapshot of the velocity field for L = 
256L/. Lower: zoom of the field in an interval of length 
10Z//, represented by dashed lines in the upper figure, (b) 
Minimizing trajectories in space time for the solution to the 
randomly forced Burgers equation periodic in space with pe- 
riod L = 256L/ and over a time interval of length T — 100. 
The various trajectories are chosen sufficiently separated at 
time t — in order to lighten the picture. The minimizers 
converge to each other in a rather nonuniform manner, giving 
rise to a tree-like structure in space-time. 
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of branches decreases and a unique trunk emerges around 
the global minimizer. Note that the same type of behav- 
ior is also observed for shocks (see Fig. 4(a) I by running 
forward in time. The intermediate time asymptotics cor- 
respond to the time scale for which there is still a large 
number of such branches. 





space 



(b) 



FIG. 4: (a) Shock trajectories for aspect ratio L/Lf =32 and 
with T = 10; the different gray areas correspond to the space- 
time domains associated to the different smooth pieces flj 's of 
the velocity field at time t = 0. (b) Sketch of the space-time 
evolution of a given smooth piece flj located between two 
shock trajectories X\(t) and Xi{t) which merge at time Tj. 
The trajectories of some fluid particles are represented; the 
shock trajectories are represented as bold lines; the dashed 
curve represents the trajectory of the unique fluid particle 
that is simultaneously absorbed by both shocks X\ and Xi 
at the time of their merger. This trajectory is a T-global 
minimizer if T > T . 



The velocity field at a given time, taken for conve- 
nience to be t = 0, consists of smooth pieces separated 
by shocks. Let us denote by {flj} the set of intervals 



in [0, L), on which the solution u(-,0) is smooth. The 
boundaries of the Qj's are the shock positions. Each of 
these shocks is associated to a root-like structure formed 
by the trajectories of the various shocks which have 
merged at times t < to form the shock under consid- 
eration (see Fig. 4(a) I. This root-like structure contains 
the whole history of the shock and in particular its age. 
Indeed, if the root has a finite depth, the shock under 
consideration has only existed for a finite time. A T- 
global shock is defined as a shock whose associated root 
is deeper than — T. One can give the dual definition for 
T-global minimizers. All the smoothness intervals Qj de- 
fined above, except that which contains the global min- 
imizer, will b e enti rely absorbed by shocks after a finite 
time (see Fig. 4(b) |. For each of these pieces, one can de- 
fine a life-time Tj as the time when the last fluid particle 
contained in this piece at time t = enters a shock. It 
corresponds to the first time for which the shock located 
on the left of this smooth interval at time t = merges 
with the shock located on the right. When the life-time 
of such an interval exceeds T, the trajectory of this latter 
fluid particle is called a T-global minimizer. 



o 
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FIG. 5: Density of T-global shocks as a function of T for three 
different system sizes L/Lf — 64, 128 and 256; average over 
100 realizations. 

Hence, for any age T, one can define a set of T-scale 
objects in the spatial domain [0, L) at a given time t. 
We define their spatial density as the number of such 
objects divided by the size of the domain L. The den- 
sity p(T) of T-global shocks is investigated numerically 
in the kicked case by using a two-step method. Firstly, 
the simulation is run up to some large time t for which 
the statistically stationary regime is reached; secondly, 
each shock present at time t is tracked backward in time 
down to the instant of its creation, giving an easy way to 
determine the density p(T). We also average numerical 
results with respect to the forcing realizations. It can be 
seen from Fig. [5] on which the density for three different 
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aspect ratios L = /Lf is presented, that the behavior of 
p(T) is independent of L, and that it displays a power- 
law p(T) oc T~ 2 / 3 in the intermediate asymptotic range 
Lf/u rms «T«L/ii 

We now present a simple phenomenological argument 
to explain the scaling exponent 2/3. Consider the solu- 
tion at a given time (t = 0, for instance). Denote by 
i(T) the typical spatial separation scale for two nearest 
T-global shocks. Obviously l(T) ~ l/p{T). The mean 
velocity of a spatial segment of length I is given by 



bi = — I u(x, 0) dx 



(3.2) 



Since the expected value (u(x, 0)) = it is natural to 
assume that for large I one has the following asymptotics 



[v,v+e] 



a (x, 0) dx ~ Vt, 



(3.3) 



which gives bt ~ l~ x l 2 for the mean velocity fluctuations. 
Consider now the rightmost minimizer corresponding to 
the left T-global shock and the leftmost minimizer related 
to the right one. Since there are no T-global shocks in 
between, it follows that the two minimizers we selected 
approach each other backward in time for times of the 
order of — T. This means that the backward-in-time dis- 
placement of a spatial segment of length 0(£) is itself 
0(1) for time intervals of the order of T. The correspond- 
ing displacement is given as the sum of two competing 
effects: the first, a kind of drift induced by the local mean 
velocity bi, is connected with the mean velocity fluctua- 
tions and gives a displacement oc b{T\ the second effect 
is a standard diffusive contribution oc T 1 / 2 related to the 
diffusive behavior of the minimizing trajectories. Taking 
both into account, we get 



i ~ dTr 1/2 + C 2 T 1/2 , 



(3.4) 



where C\ and C 2 are numerical constants. It is easy to 
see that the dominant contribution comes from the first 
term. Indeed, if the second term dominates, i.e. I oc T 1 / 2 , 
the first term then gives a contribution of the order of 
£ 3 / 2 which contradicts to l|3.4|l . Hence, one should have 
I ~ CiT£ -1 / 2 , leading to the scaling behavior 



£(T) oc T 2/3 , p(T) oc T~ 2/3 . 



(3.5) 



As we have already discussed above, T-global shocks are 
shocks older than T. Denote by p(A) the PDF for the 
age of shocks. More precisely, p(A) is a density in the 
stationary regime of a probability distribution of the age 
A{t) of a shock, say the nearest to the origin. It follows 
from H3.5f) that the probability of shocks whose age is 
larger than A decays like A~ 2 / 3 , which implies the fol- 
lowing asymptotics for the PDF p(A): 



IV. TWO RELATED STATISTICAL MODELS 

We next discuss a very close analogy between the problem 
just considered and the dynamical behavior of a simple 
first-order stochastic ODE. Consider the following one- 
dimensional stochastic Ito equation: 



dX = f(X)dW(t), 



(4.1) 



where < C\ < f(x) < C% < oo, x £ M. We are in- 
terested in the stochastic flow generated by 1)4. In 
other words, the aim is to understand the dynamical 
properties in the large-time asymptotics of the trajec- 
tories corresponding to the solutions of (|4.1|) for different 
initial values. Let Xi(t) and X2(t) with t > be the 
two solutions associated to the initial values Xi(0) = xi 
and X 2 (0) — x 2 . It is well-known that the difference 
[Xi (t) — X2 (t)) is a martingale which does not change its 
sign; this implies that the limit 



c( Xl ,x 2 )= lim [X x {t) - X 2 (t)] 

t — >oo 



(4.2) 



p(A) oc A -5 / 3 . 



(3.6) 



exists and is finite. 

Consider first the case when f(x) is periodic. It is easy 
to see that c(x\,x 2 ) has to be equal to a period of the 
function f(x) (see, e.g., Fig. 6(a) I. Let, for simplicity, the 
smallest period of f(x) be equal to 1. Then, c(xi,x 2 ) 
can take only two values: either or 1. In fact, it is 
possible to show (see 0,EEH3) that the following 
picture holds. There exists a random periodic sequence 
of points x(k) — x*(uj) + fc, k € Z, x* G [0, 1) such that 
c(xi,x 2 ) = for all £1,2:2 € + k, x*(uj) + k + 1) 

and c(xi,x 2 ) = 1 for x\ € (x* (u) + k — 1, x*(u) + fc) and 
x 2 e (x^(cu) + k, x*(uj) + k+1). Moreover, in both cases, 
the convergence is exponential: 

\c(x 1 ,x 2 )-(X l (t)-X 2 (t))\ ~Cexp(-At), (4.3) 

where A > is a non-random Lyapunov exponent (see 
Fig. 6(b) I. The random variable x*(w) which determines 
the positions of the separation points x(k) depends on 
the realization ui of the white noise dW(t). If we con- 
sider the stochastic flow on the unit circle, then all the 
trajectories are stable except the one originating from 
x*(oj). The stability is governed by the non-random pos- 
itive Lyapunov exponent A. In contrast, the trajectory 
which starts at x*(uj) is unstable and pushes nearby tra- 
jectories away. Hence, x*(uj) plays a role similar to that 
of the global minimizer for the Burgers equation with 
a spatially periodic force. An exponentially small neigh- 
borhood of it gets mapped into the whole unit circle apart 
from another exponentially small interval which is the 
image of the rest of the circle. This instability island 
around x*(lj) is similar to the interval of fluid particles 
which are not absorbed by shocks until time t. 

All the results above in the case of periodic f(x) have 
been well-known for about 20 years and were presented 
here just to emphasize an analogy with the main shocks 
and minimizers for Burgers equation in a periodic sit- 
uation. At the same time, the dynamical behavior in 
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FIG. 6: Behavior of the solutions to the stochastic ODE 14.111 
for the periodic function f(x) — 2 + sin a;, (a) Unwrapped 
trajectories of the solutions converging to one of the periodic 
image, (b) Time behavior, for a given realization, of the sep- 
aration between two trajectories xi(t) and X2(t) associated to 
two different initial conditions. The two trajectories approach 
each other exponentially fast with a non-random rate A. 

the non-periodic case is more complicated and much less 
studied. As we explain below, this behavior is also very 
similar to the dynamical picture governed by the re- 
pelling T-global shocks in the case of non-periodic Burg- 
ers equation. Since f(x) is non-periodic c(x%, x 2 ) = and 
all trajectories are asymptotic to each other as t — > oo. 
However, as in the case of non-periodic Burgers equation, 



this convergence is very non-uniform. Although there are 
no globally unstable positions like x(k) — x*{uj) + k, for a 
given time scale T ^> 1, one still has exponentially small 
T-instability islands which separate different intervals of 
stability from each other. For a given time T, the image 
under the stochastic flow of the T-instability islands cov- 
ers almost all the real line apart from exponentially small 
pieces which are images of the stability intervals. Those 
exponentially small pieces, which we call T '-stability in- 
tervals play a role similar to that of the T-global shocks 
in Burgers turbulence. For a larger time scale Ti > T a 
majority of T-instability island will become stable, while 
some of them will still be unstable. Those islands are 
in fact Ti-unstable. One can ask the same question as 
studied above about the asymptotics of the density of 
T-instability islands in the limit as T — > oo. Our pre- 
liminary numerical results indicate that the density of 
T-instability islands scales as T -1 / 2 . Notice that this 
scaling is different from the scaling T -2 / 3 in the case 
of Burgers equation. We believe that in the case of a 
stochastic flow the T-instability islands and the density 
of them can be defined and studied rigorously. The sim- 
plification is connected with the absence of a slow drift 
which one has for Burgers equation due to fluctuations 
of the average velocity. In the case of the stochastic flow 
the whole asymptotic dynamical picture is formed by just 
two factors: diffusion and hyperbolic contraction. As a 
result, the rigorous analysis of the model seems to be a 
realistic task. 

Wc finally suggest a simple statistical model which cap- 
tures the essential characteristics of the dynamics and the 
interactions of T-global shocks and T-stability intervals. 
We believe that this model is also quite interesting on 
its own. Consider an infinite system of particles on a 
one-dimensional integer lattice Z. Two different parti- 
cles cannot occupy the same site, and some sites are in 
general not occupied. Each particle is associated to an 
integer age A > 0. The model consists in describing a dis- 
crete time evolution of particles. To get a configuration 
of particles at the next step, we proceed as follows. First 
of all, with a probability < p < 1, we generate indepen- 
dently new particles in every non-occupied site and we 
set the age A = for all newly born particles. Then, each 
particle independently jumps, with a probability 1/2, ei- 
ther to the right or to the left by a distance 1/2 and we 
add 1 to the age of all particles. Notice that, after one 
time step, the particles are located on the sites of the dual 
lattice. If there are two particles which jumped into the 
same point of the dual lattice, we then assume that the 
older one absorbs the younger, that is we ascribe to this 
site a particle whose age is the maximum of the two ages 
of the merging particles. It is easy to sec that this dynam- 
ics reaches a statistical steady state and one can define 
an integer-valued random variable A which is the age of 
the particle closest to the origin, say from the right. The 
probability distribution of this random variable, namely 
the asymptotics as A — > oo of the probability p{A), is 
closely connected with the statistics of T-instability in- 
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tervals for the non-periodic stochastic flow which we dis- 
cussed above. A numerical study of the model suggests 
that p(A) scales like A~ 3 / 2 which corresponds to the scal- 
ing T^ 1 ' 2 for the density of T-instability islands. We 
postpone the detailed numerical and theoretical analy- 
sis of the non-periodic stochastic flow and the statistical 
model of interacting particles until a future publication. 



V. CONCLUDING REMARKS 

We have studied the inviscid randomly forced Burgers 
equation with non-periodic forcing on the whole real line 
started at t = — oo. Our results indicate the existence of 
stationary regime which corresponds to the velocities of 
one-sided minimizers and suggest the following picture. 
At any given time t (say t — 0) and any given x € K, 
there exists at least one one-sided minimizer. However, 
due to fluctuations of the positions of one-sided mini- 
mizers, there are no global minimizers. This means that 
any fluid particle gets absorbed by a shock after a cer- 
tain time. Any two one-sided minimizers are asymptotic 
to each other backwards-in-time, that is the distance be- 
tween them tends to zero as t tends to — oo. However, 
this convergence is very nonuniform. For fixed x, y G R 
one-sided minimizers which originated at t — in x and 
y will approach each other exponentially fast beyond a 
certain correlation time T(x,y). This correlation time 
is of the order of the maximum T for which a T-global 
shock exists between x and y. One can say that T-global 
shocks form a hierarchical structure separating different 
stability domains from each other. By stability domains 
we understand intervals for which one-sided minimizers 
converge exponentially fast with a correlation time of the 
order of the turnover time Lf/u rms . The separation by 
T-global shocks form separation "walls" between these 
different stability intervals for times of the order of — T. 
For larger backward times, one-sided minimizers from the 
neighboring stability intervals are exponentially asymp- 
totic to each other. Another interpretation of T-global 
shocks is connected with their time of creation. Every 
shock can be traced backward only for a finite time in- 
terval. For T-global shocks, this time interval which de- 
termines the age of the shock is larger than T. Since 
all shocks have a finite age, it follows that there are no 
true main shocks in the non-periodic situation. Our re- 
sults suggest that the large-T asymptotics of the density 
of T-global shocks follows the power-law p(T) oc T~ 2 / 3 . 
This gives a power-law behavior with exponent 5/3 for 
the PDF of the age of shocks: p(A) oc A -5 / 3 . Another 
related exponent is connected with the absorption times, 
defined in the following way. The absorption time T(£) 
is the time after which all fluid particles in the spatial 
interval [—£/2, £/2] are absorbed by shocks. Using the 
duality between the forward-in-time behavior of fluid par- 
ticles (minimizers) and the backward-in-time behavior of 
shocks, we find that the absorption times T(l) oc £ 3 / 2 at 
large i. 




10 1 io° io' io 2 io 3 io 4 io 5 



wavenumber 271 k 

FIG. 7: Spectrum (^(k)) of the velocity potential in the 
stationary regime for the aspect ratio L/Lf = 128. This 
spectrum contains two power-law ranges: at wavenumbers 
k 2> L/Lf, the traditional oc k~ 4 inertial range connected to 
the presence of shocks in the solution and, for k <C L/Lf, an 
"inverse cascade" oc k~ 2 associated to the large-scale fluctu- 
ations of tp 

Notice that the power-law behavior of the density p(T) 
of T-global shocks can be interpreted in term of a kind of 
inverse cascade in the spectrum of the solution (although 
there is no conserved energy-like quantity). Indeed, the 
fluctuations (|3.3fl of the mean velocity suggest that, for 
large-enough separations i, the velocity potential incre- 
ment scales like 

\ip(x + £, t)-ip(x, t)\ ocf 1 / 2 . (5.1) 

This behavior is responsible for the presence of an inter- 
mediate power-law range with exponent —2 in the spec- 
trum of the velocity potential at wavenumbers smaller 
than the forcing scale (see Fig.[7|). This power-law range 
corresponds for the velocity spectrum to an equipartition 
of kinetic energy, meaning that the large-domain asymp- 
totics may be considered in the same universality class as 
the KPZ problem for interface dynamics. It is also im- 
portant to notice that, in order to obtain this —2 range at 
small wavenumbers, the spectrum of the forcing poten- 
tial must go to zero faster than k~ 2 as k — + 0. Otherwise, 
the leading behavior is non-universal and depends on the 
functional form of the forcing correlation. 

The randomly forced Burgers equation in an un- 
bounded domain with a different type of forcing was also 
considered in Ref. [2l|, where it was assumed that the 
forcing potential has at any time its global maximum 
and its global minimum in a prescribed compact part of 
the space. Such a forcing leads to a completely different 
behavior. In particular, in this case there exists a global 
minimizer located in a finite spatial interval for all times 
and all other minimizers are asymptotic to it in the limit 
t — » oo. 
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We finally stress that, in this paper, we have mostly 
discussed the one-dimensional case. It is of particular 
interest to analyze the effects of non-compactness of the 
domain and the intermediate time asymptotics in higher 
dimensions. The notion of main shock is then replaced 
by that of topological shock 0, ^} , which are no more 
isolated points, but spatially extended objects. The nat- 
ural problem in this setting is to study geometrical and 
statistical properties of the T-global shocks and to find 
the asymptotic behavior of their age distribution. 
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